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(Dated: February 7, 2008) 

We present measurements of the A° lifetime in the exclusive decay channel A^ — > J/ipA , with 
J/tp —* an d A — > p-K~ , the B° lifetime in the decay B° — > J/tpKg with J/ip — > fi + /i~ 

and Kg — > tv + tt~, and the ratio of these lifetimes. The analysis is based on approximately 250 
pb _1 of data recorded with the D0 detector in pp collisions at -^3=1.96 TeV. The A° lifetime is 
determined to be r(A°) = 1.22±°; 22 .(stat) ± 0.04(syst) ps, the B° lifetime t(B°) = 1.40±°, ^(stat) ± 
0.03(syst) ps, and the ratio t(A°)/t(B°) = O.STt^JKstat) ± 0.03(syst). In contrast with previous 
measurements using semileptonic decays, this is the first determination of the A° lifetime based on 
a fully reconstructed decay channel. 



PACS numbers: 14.20.Mr, 14.40.Nd, 13.30.Eg, 13.25.Hw 



Calculations based on a simple quark-spectator 
model predict that the lifetimes of all b hadrons are 
equal. When non-spectator effects are taken into ac- 
count, they give rise to a lifetime hierarchy of t(B + ) > 



t{B°) « t(B° s ) > r(A°) > t(B+) '% Measurements 
of 6-hadron lifetimes therefore provide means to deter- 
mine the importance of non-spectator contributions in 
6-hadron decays. For comparison with theory, measure- 
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merits of lifetime ratios are preferred over individual life- 
times. The ratio t(A^)/t(B°) has been the source of 
theoretical study since early calculations [j| predicted a 
value greater than 0.90, almost two sigma away from the 
measurement average |4j, 0.800 ± 0.053. Recent calcu- 
lations Q of this ratio that include higher order effects 
have reduced this difference. A current compilation 0| 
of lifetime ratio data for b hadrons is consistent with 
early calculations Q that include non-spectator effects. 
Previous measurements of r(A°) used semileptonic decay 
channels that suffer from uncertainties arising from un- 
detected neutrinos. A measurement of the lifetime using 
fully reconstructed A^ decays is free from the neutrino 
ambiguities. The Tevatron Collider at Fermilab is the 
only operating accelerator where A° baryons are being 
produced and studied. 

In this Letter we report a measurement of the A° life- 
time in the decay channel A^ — > J/ipA Q , and its ratio to 
the B° lifetime from the B° — > J/ipK s decay channel. 
This B° decay channel is chosen because of its similar 
topology to the A° decay. The J/%p is reconstructed in the 
decay mode, the A in pn~ , and the Kg in tt + it~; 
throughout this Letter the appearence of a specific charge 
state will also imply its charge conjugate. The data used 
in this analysis were collected during 2002-2004 with 
the D0 detector in Run II of the Tevatron Collider at 
a center-of-mass energy of 1.96 TeV, and correspond to 
an integrated luminosity of approximately 250 pb _1 . 

The components of the D0 detector most relevant 
for this measurement are the charged-particle tracking 
systems and the muon detector. The D0 tracker consists 
of a silicon microstrip tracker (SMT) and a central fiber 
tracker (CFT) that are surrounded by a superconducting 
solenoid magnet that produces a 2 T central magnetic 
field. The SMT has approximately 800,000 individual 
strips, with a typical pitch of 50 — 80 /im, and a design op- 
timized for tracking and vertexing capability for |?7| < 3 
(r/ = — ln[tan(6'/2)] and is the polar angle). The sys- 
tem has a six-barrel longitudinal structure interspersed 
with sixteen disks. Each barrel consists of four layers ar- 
ranged axially around the beam and the disks are placed 
perpendicular to the beam. The CFT has eight thin coax- 
ial barrels, each supporting two doublets of overlapping 
scintillating fibers of 0.835 mm diameter, one doublet 
being parallel to the collision axis, and the other alter- 
nating by ±3° relative to the axis. For charged particles, 
the resolution for the distance of closest approach as pro- 
vided by the tracking system is approximately 50 /im for 
tracks with px ~ 1 GeV/c, and improves asymptotically 
to 15 /im for tracks with px > 10 GeV/c, px is the compo- 
nent of the momentum perpendicular to the beam pipe. 
Preshower detectors and electromagnetic and hadronic 
calorimeters surround the tracker. A muon system is lo- 
cated beyond the calorimeter, and consists of multilayer 
drift chambers and scintillation trigger counters preced- 
ing 1.8 T toroidal magnets, followed by two similar layers 



beyond the toroids. Muon identification for |?7| < 1 relies 
on 10 cm wide drift tubes, while 1 cm mini-drift tubes 
are used for 1 < |?/| < 2. 

Primary vertex (PV) candidates are determined for 
each event by minimizing a \ 2 function that depends 
on all the tracks in the event and a term that represents 
the beam spot constraint. The beam spot is the run- by- 
run average beam position, where a run typically lasts 
several hours. The beam spot is stable during the peri- 
ods of time when the proton and antiproton beams are 
kept colliding continuously and can be used as a con- 
straint for the primary vertex fit. The initial primary 
vertex candidate and its \ 2 are obtained using all tracks. 
Next, each track used in the \ 2 calculation is removed 
temporarily and the x 2 is calculated again; if the \ 2 de- 
creases by 9 or more, this track is discarded from the PV 
fit. This procedure is repeated until no more tracks can 
be discarded. Additional primary vertices are obtained 
by applying the same algorithm to the discarded tracks 
until no more vertices are found. 

We base our data selection on defined objects such 
as charged tracks and muons. Although we do not re- 
quire any specific trigger to select our sample, most of 
the events selected fire dimuon or single muon triggers. 
Preliminary selection of dimuon events requires the pres- 
ence of at least two muons of opposite charge recon- 
structed in the tracker and the muon system. We re- 
quire that at least one of the muon candidates consists 
of a central track, defined by hits in the SMT and CFT, 
matched with muon track segments in all three layers of 
the muon system. For the second muon, we require a 
central track matched with hits in at least the innermost 
layer of the muon system. The sample of J /ip — * 
candidates consists of events with at least two muons, 
with trajectories constrained in a fit to a common vertex. 
The fit must have a % 2 probability greater than 1%, and 
the invariant mass of the dimuons must be in the range 
2.80 < Mfifi < 3.35 GeV/c 2 . To reconstruct A° and B° 
candidates, the J/tp events are examined for A and K s 
candidates. The A — > pir~ candidates are required to 
have two tracks of opposite charge which must originate 
from a common vertex with a \ 2 probability greater than 
1%. A candidate is selected if the mass of the proton- 
pion system after the vertex-constrained fit falls in the 
1.100 < M p7r < 1.128 GeV/c 2 window. The proton mass 
is assigned to the track of higher momentum, and the px 
of the A is required to be greater than 2.4 GeV/c. The 
K s — > 7r + 7r~ selection follows the same criteria, except 
that the mass window is 0.460 < M„ < 0.525 GeV/c 2 , 
andpT >l-8 GeV/c. 

We reconstruct the A° and B° by performing a con- 
strained fit to a common vertex for either the A or K s 
and the two muon tracks, with the latter constrained to 
the J/tp mass of 3.097 GeV/c 2 Q. Because of their long 
decay lengths, a significant fraction of A and K s will 
decay outside the SMT. Therefore, to maintain good ef- 
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ficiency, no SMT hits are required on the tracks of the 
decay particles. To reconstruct the A° (B°), we first find 
the A (K s ) decay vertex, and then extrapolate the mo- 
mentum vector of the ensuing particle and form a vertex 
with it and the two muon tracks belonging to the J/ip. 
The precision of the A° (B°) vertex position is dominated 
by the two muon tracks from the J/tp- If more than one 
candidate is found in the event, the candidate with the 
best x 2 probability is selected as the A° (B°) candidate. 

We determine the lifetime of a A° or B by measur- 
ing the distance traveled by each 6-hadron candidate in a 
plane transverse to the beam direction, and then applying 
a correction for the Lorentz boost. We define the trans- 
verse decay length as L xy — L xy ■ Pt/pt where L xy is 
the vector that points from the primary to the secondary 
vertex and pt is the transverse momentum vector of the 
b hadron. The proper decay length (PDL) for a 6-hadron 
candidate is then given by: 

P r>T - Lxy - T cMb h ) 

where (j3^f) B , and Mb are the transverse boost and the 
mass of the b hadron, respectively. In our measurement, 
the value of Mr in Eq. Q]is set to the PDG mass value 
of Ag or B° In our final selection of Ag and B° 
candidates, we require an error of less than 100 /im on 
the PDL and we also require a total momentum greater 
than 5 GeV/c. 

We perform an unbinned likelihood fit to measure the 
A° and B° lifetimes. The inputs for the fit are the mass, 
PDL and PDL error of the candidates. Candidates with 
invariant masses in the range of 5.1 to 6.1 GeV/c 2 for the 
Ag and 4.9 to 5.7 GeV/c 2 for the B° are selected; these 
ranges include sideband regions that are used to model 
the PDL distributions of backgrounds. The likelihood 
function, £, is defined by: 

£ = U.f=i[fsSM(M j )S L (Xj,a j ) {) 
+ (1- fs)B M (M J )B L (X J ,a J )}, V ' 

where N is the total number of selected events, f s is 
the fraction of signal events in the sample, Sm and Bm 
are the probability distribution functions used to model 
the mass distributions for signal and background, respec- 
tively, and Sl and Bl model the distributions of proper 
decay lengths for signal and background. The mass for 
signal is modeled by a Gaussian distribution and the mass 
for background is described by a second-order polyno- 
mial. The PDL distribution for signal is described by the 
convolution of an exponential decay, whose decay con- 
stant is one of the parameters of the fit, with a resolution 
function represented by a single Gaussian function: 

G(A -^ ) = ^ eXP (^)' (3) 




Invariant Mass (GeV/c 2 ) 



FIG. 1: Invariant mass distribution for Aj; candidate events. 
The points represent the data, and the curve represents the 
result of the fit. The mass distribution for the signal is shown 
in gray. 



where Aj and Oj represent the PDL and its error respec- 
tively for a given event j, and the s parameter is intro- 
duced in the fit to account for a possible misestimate of 
cry. The PDL distribution for background is described 
by a sum of a resolution function representing the zero- 
lifetime component, negative and positive exponential 
decay functions modeling combinatorial background, and 
an exponential decay that accounts for long-lived heavy 
flavor decays. We minimize —2 ln£ to extract the param- 
eters: ct(A°) = 366±g a™ and ct(B°) = 419±H fxm. 
From the fits, we get s = 1.27 ±0.10 and s = 1.39 ±0.05 
for the A° and B° respectively; and the number of signal 
events 61 ± 12 A° and 291 ± 23 B°. Figs [T] and (Figs 
and 15 show the mass and proper decay length distribu- 
tions for the A° (B°) candidates, respectively, with the 
results of the fits superimposed. 

Table[l]summarizes the systematic uncertainties on our 
measurements. The contribution from the uncertainty on 
the detector alignment is estimated by reconstructing the 
B° sample with the positions of the SMT sensors shifted 
outwards radially by the alignment error in the radial 
position of the sensors and then fitting for the lifetime. 
We estimate the systematic uncertainty due to the reso- 
lution on the PDL by using two Gaussian functions for 
the resolution model. The contribution to the system- 
atic uncertainty from the model describing background 
PDLs is studied by varying the parametrizations of the 
different components: (i) the exponential functions are 
replaced by exponentials convoluted with the resolution 
function of Eq. |3J (ii) a uniform background is added 
to account for outlier events (this has only a negligible 
effect), and (iii) the positive and negative short-lived life- 
time components are forced to be symmetric. To study 
the systematic uncertainty due to the model for the mass 
distributions, we vary the shapes of the mass distribu- 
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FIG. 2: Distribution of proper decay length for candi- 
dates. The points are the data, and the solid curve is the sum 
of the contributions from signal (gray) and the background 
(dashed-dotted line). 




-0.1 -0.05 



0.05 0.1 0.15 0.2 0.25 0.3 
Proper decay length (cm) 



FIG. 4: Distribution of proper decay length for B° candi- 
dates. The points are the data, and the solid curve is the sum 
of the contributions from signal (gray) and the background 
(dashed-dotted line). 
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FIG. 3: Invariant mass distribution for B° candidate events. 
The points represent the data, and the curve represents the 
result of the fit. The mass distribution for the signal is shown 
in gray. 



tions for signal and background. For the signal, we use 
two Gaussian functions instead of a single one, and for 
the background distribution, a linear function instead of 
the nominal quadratic form. 

The lifetime of the long-lived component of the back- 
ground varies with mass. This results in an uncertainty 
in the decay constant of the background under the mass 
peaks. We obtain the systematic uncertainty due to this 
effect by modelling the long-lived background with two 
exponentials instead of a single exponential. The decay 
constant of one of the two exponentials is determined 
from a fit in the low-mass sideband, and the other de- 
cay constant is determined from the high-mass sideband. 
The low-mass sideband is defined as the mass window 
4.900-5.149 GeV/c 2 for B° and 5.100-5.456 GeV/c 2 for 



Ag and the high-mass sideband as 5.389-5.700 GeV/c 2 
and 5.768-6.100 GeV/c 2 respectively. We perform the 
fit incorporating the linear combination of exponentials 
with the decay constants fixed to the values obtained in 
the low- and high-mass sidebands fits and allowing the 
coefficients of the linear combination to float. The sys- 
tematic uncertainty quoted is the difference between the 
values we get from this fit and the nominal. 

We also study the contamination of the A° sample by 
B° events that pass the A° selection. From Monte Carlo 
studies, we estimate that 19 B° events are reconstructed 
as A° events. The invariant masses of the B° events en- 
tering the A° sample are distributed almost uniformly 
across the entire mass range, and do not peak at the 
A° mass. Their proper decay lengths therefore tend to 
be incorporated in our model of the long-lived heavy- 
flavor component of the background. To estimate the 
systematic uncertainty due to this contamination, we fit 
the mass and proper decay length distributions of the 
misidentified events in the MC samples, add this con- 
tribution to the likelihood with fixed parameters, and 
perform the fit again. The difference between the two re- 
sults is quoted as the systematic uncertainty due to the 
contamination. 

The fitting procedure is tested for the presence of bi- 
ases by generating 1000 Monte Carlo experiments, each 
with the same statistics as our data samples. For the 
generated events, the PDL errors are taken from data, 
and the mass and PDL distributions are described by 
the probability distribution functions used in data, with 
parameters obtained from the fit. The fits performed on 
these Monte Carlo experiments indicate that there is no 
bias inherent in the procedure. 

We also perform several cross-checks of the lifetime 
measurements. In particular, a fit is done where the back- 



7 



TABLE I: Summary of systematic uncertainties in the mea- 
surement of ct for A° and B° and their ratio. The total un- 
certainties are also given combining individual uncertainties 
in quadrature. 



Source 




B° (>m) 


Ratio 


Alignment 


5.4 


5.4 


0.002 


Model for PDL resolution 


6.7 


2.7 


0.010 


Model for PDL background 


2.7 


3.1 


0.005 


Model for signal mass 


0.2 


0.0 


0.000 


Model for background mass 


2.5 


6.2 


0.007 


Long-lived components 


1.5 


0.1 


0.003 


Contamination 


8.8 


0.8 


0.023 


Total 


12.9 


9.2 


0.028 



ground is modeled using only sideband regions, the J /ip 
vertex is used instead of the 6-hadron vertex, the mass 
windows are varied, the reconstructed &-hadron mass is 
used instead of the Particle Data Group jj] value, and 
the sample is split into different pseudorapidity regions 
or different regions of azimuth. All results obtained with 
these variations are consistent with our central values. 

The results of our measurement of the A° and B° life- 
times are summarized as: 

r(Ag) = 1.22±g;?| (stat) ± 0.04 (syst) ps, (4) 
t(B°) = lAOto.ll (stat) ±0.03 (syst) ps. 

These can be combined to determine the ratio of life- 
times: 



= 0.87±°;£ (stat) ± 0.03 (syst), (5) 

where we determine the systematic uncertainty of the 
ratio by varying each parameter in the two samples si- 
multaneously and quoting the deviation in the ratio as 
the systematic uncertainty due to that source. 

In conclusion, we have measured the A° lifetime in the 
fully reconstructed exclusive decay channel J /iphP . This 
is the first time that this lifetime has been measured in an 
exclusive channel. The measurement is consistent with 
the world average, 1.229 ± 0.080 ps 4.], and the Ag to 
B° ratio of lifetimes is also consistent with theoretical 
predictions 0. 0. Q 
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